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Abstract 

A formula for the number of toroidal m x n binary arrays, allowing 
rotation of the rows and/or the columns but not reflection, is known. Here 
we find a formula for the number of toroidal mxn binary arrays, allowing 
rotation and/or reflection of the rows and/or the columns. 
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1 Introduction 



The number of necklaces with n beads of two colors when turning over is not 
allowed is 

i5>(d)2" /d , (1) 



n 
d 



where tp is Euler's phi function. When turning over is allowed, the number 
becomes 

^£^)^ + ( 2( "T 2 ifnisodd ' (2) 

2n V I 3 • 2"/ 2 - 2 if n is even. V 7 

d\n \ 

These are the core sequences A000031 and A000029, respectively, in [2]. 

Our concern here is with two-dimensional versions of these formulas. We 
consider an m x n binary array. When opposite edges are identified, it becomes 
what we will call a toroidal binary array. Just as we can rotate a necklace 
without effect, we can rotate the rows and/or the columns of such an array 
without effect. The number of (distinct) toroidal m x n binary arrays is 

E^ c M d ) 2m " /lcmM) > ( 3 ) 

c | m d | n 

where 1cm stands for least common multiple. This is A184271 in [2]. The 
diagonal is A179043. Rows (or columns) 2-8 are A184264-A184270. Row (or 
column) 1 is of course A000031. 

Our aim here is to find the formula that is related to (3) in the same way that 

(2) is related to (1). More precisely, we wish to count the number of toroidal 
m x n binary arrays allowing rotation and/or reflection of the rows and/or the 
columns. At present, the resulting sequence, the diagonal, and the rows (or 
columns) other than the first one, are not found in [2] . Row (or column) 1 is of 
course A000029. 

For an alternative description, we could define a group action on the set of 
m x n binary arrays, which has 2 mn elements. If the group is C m x C n , where 
C m denotes the cyclic group of order m, then the number of orbits is given by 

(3) (see Theorem 1 below). If the group is D m x D n , where D m denotes the 
dihedral group of order 2m, then the number of orbits is given by Theorem 2 
below. 

To help clarify the distinction between the two group actions, we provide an 
example. There is no distinction in the 2x2 case, so we consider the 3x3 case. 
When the group is C3 x C3 (allowing rotation of the rows and/or the columns 
but not reflection), there are 64 orbits, as shown in Table 1. 

When the group is D3 x D3 (allowing rotation and/or reflection of the rows 
and/or the columns), there are 36 orbits, as shown in Table 2. 

Our interest in the number of toroidal m x n binary arrays allowing rotation 
and/or reflection of the rows and/or the columns derives from the fact that this 
is the size of the state space of the projection of the Markov chain in [1] under 
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Table 1 : A list of the 64 orbits of the group action given by the group C3 x C3 
acting on the set of 3 x 3 binary arrays. (Rows and/or columns can be rotated 
but not reflected.) Each orbit is represented by its minimal element in 9-bit 
binary form. Bars separate different numbers of Is. 
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Table 2: A list of the 36 orbits of the group action given by the group D 3 x D 3 
acting on the set of 3 x 3 binary arrays. (Rows and/or columns can be rotated 
and/or reflected.) Each orbit is represented by its minimal element in 9-bit 
binary form. Bars separate different numbers of Is. 
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the mapping that takes a state to the orbit containing it. This reduction of the 
state space, from 512 states to 36 states in the 3x3 case for example, makes 
computation easier. 



2 Rotations of rows and columns 

Let X m , n := {0, i}{o,i,-,m-i}x{o,i,...,n-i} bc the set f m x n arrays f o s and 
Is, of which there are 2 mn . Let a(m, n) denote the number of orbits of the group 
action on X m>n by the group C m x C n . In other words, a(m, n) is the number 
of (distinct) toroidal m x n binary arrays, allowing rotation of the rows and/or 
the columns but not reflection. 

Theorem 1. 

a(m, n) = ^ E E ^( c M d ) 2" m / lcm ( c < d > . (4) 

c | 77i d | n 

Proof. By the Polya enumeration theorem [3], 

1 rn— 1 n— 1 

^^-EE 2Ay . (5) 

i=0 j=0 

where Aij is the number of cycles in the permutation crV- 7 ; here a rotates the 
rows (row becomes row 1, row 1 becomes row 2, . . . , row m — 1 becomes row 
0) and t rotates the columns. For example, Aqq = mn because the identity 
permutation has mn fixed points, each of which is a cycle of length 1. 

It is well known that, if d divides n, then the number of elements of C n that 
are of order d is tp(d). So if c divides m and d divides n, then the number of 
pairs (i, j) such that a 1 is of order c and t j is of order d is <p(c)ip(d). And if a 1 is 
of order c and t 3 is of order d, then <j 1 t 3 is of order lcm(c, d) because a 1 and t 3 
commute. Consequently, lcm(c, d) is the length of each cycle of the permutation 
a % T J , so = mn/\cm(c,d), and (4) follows from (5). □ 

Clearly, a(l, n) reduces to (1). Table 3 provides numerical values of a(m, n) 
for small m and n. 



Table 3: The number a(m, n) of toroidal m x n binary arrays, allowing rotation 
of the rows and/or the columns but not reflection, for m,n = 1, 2, . . . , 8. 



14 

• 10 
108 
362 
1182 
4150 



14 

64 
352 
2192 
14624 
99880 
699252 



40 
352 
4156 
52488 
699600 
9587580 
134223976 



108 
2192 
52488 
1342208 
35792568 
981706832 
27487816992 



1 1 
362 
14624 
699600 
35792568 
1908897152 
104715443852 
5864063066500 



20 
1182 
99880 
9587580 
981706832 
104715443852 
11488774559744 
1286742755471400 



36 
4150 
699252 
134223976 
27487816992 
5864063066500 
1286742755471400 
288230376353050816 
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3 Rotations and reflections of rows and columns 



Let 6(to, n) denote the number of orbits of the group action on X m ^ n by the 
group D m x D n . In other words, 6(m, n) is the number of (distinct) toroidal 
m x n binary arrays, allowing rotation and/or reflection of the rows and/or the 
columns. 

Theorem 2. 

6(to, n) = b\ (to, n) + 62(771, n) + 63(771, n) + 64(771, n), (6) 

where 



bl{m ,n) = -L- 5>(cMd)2' 



^mn/lcm(c,d) 

c I m d\n 



b 2 (m, n) 

= r (4n) -l 2 (m+l)n/2 7/7777^0' 1 ^^mn/d 

I (8n)"M2 mn/2 + 2(" i+2 W 2 l if mi s even 4n f-J, 

J(4n)- 1 ^'[2( m + 1 )Scd(j,n)/2 _ 2 mgcd(3,»»)] j/ m j s 0£ fcZ 

+ {(871)- 1 ^'[2 m s cd (^«)/ 2 + 2(™+ 2 )g cd (^«)/ 2 - 2 r ™s cd (^") +1 ] i/m is even 

With = El<j<n-l:n/gcdO» is odd' 

6 3 (to,7i) = 62(77,777), 

£7776? 

{2(mn-3)/2 ^ m flnc / n are ^ 
3 • 2 m "/ 2 ~ 3 if m and 77 raaiie opposite parity, 
7 • 2 m "/ 2-4 if 777 and 77 are even. 



Proof. Again by the Polya enumeration theorem [3] , 

m — 1 n — 1 

4mn 



m— 1 n— 1 

6(777, 77) = ^ ^ [2 A - + 2 s - + 2 C - + 2 D - ] , (7) 

»=o j=0 

where (rcsp., Bij, C,j, Z),j) is the number of cycles in the permutation erV^ 
(resp., <3 % t3 p, a l T^d, a % r^ p9); here a rotates the rows (row becomes row 1, 
row 1 becomes row 2, . . . , row 777 — 1 becomes row 0), r rotates the columns, 
p reflects the rows (rows and m — 1 are interchanged, rows 1 and to — 2 are 
interchanged, rows L m /2J — 1 and to — [ro/2j are interchanged), and 9 
reflects the columns. 

By the proof of Theorem 1, we know the form of Aij, and this gives the 
formula for 61 (to, n). 

Next we find (Bio), the entries in the 0th column of matrix B. For i = 
0, 1, . . . , m — I, the permutation a l p can be described by its effect on the rows 
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of {0, 1, . . . , m — 1} x {0, 1, . . . , n — 1}. It reverses the first m — i rows and 
reverses the last i rows. Since the reversal of k consecutive integers has k/2 
transpositions if k is even and (k — l)/2 transpositions and one fixed point if 
k is odd, the permutation of {0, 1, . . . ,m — 1} induced by a l p has (m — l)/2 
transpositions and one fixed point if m is odd, and m/2 transpositions if i is 
even and m is even, and (m — 2)/2 transpositions and two fixed points if i is odd 
and m is even. These numbers must be multiplied by n for the permutation a 1 p 
of {0, 1, ... ,m — 1} x {0, 1, ... ,n — 1}. The results are that B M = (m + l)n/2 
if m is odd, B,o = mn/2 if i is even and m is even, and B^ = (m + 2)n/2 if i 
is odd and m is even. Therefore, (4mn)^ 1 YJ!=o 2S, ° = (4n)- 1 2( m + 1 )™/ 2 if m 
is odd, whereas (4mn)" 1 X^o" 1 2S *° = (Sn)- 1 ^" 1 "/ 2 + 2( m+2 W 2 ] if m is even, 
and this gives the first term in the formula for 62 (m, n). 

We turn to Bij for i = 0, 1, . . . , m — 1 and j = 1, 2, . . . , n — 1. First, by 
a property of cyclic groups, has order d := ?i/gcd(j, n). If d is even, then, 
since a 1 p has order 2 (see the preceding paragraph), a % T 3 p has order d and all 
of its cycles have length d. In this case, B^ = mn/d = mgcd(j,n). Suppose 
then that d is odd. There are three cases: (i) m odd, (ii) i even and m even, 
and (m) i odd and m even. Recall that a l p reverses the first m — i rows 
and reverses the last i rows. In case (i), one row is fixed by <r l p, so cycles of 
a l r^ p in this row have length d and all others have length 2d. We find that 
B^ = n/d + (m - l)n/(2d) = (m + i)n/(2d) = (m + 1) gcd(j, n)/2. In case 
(ii), no rows are fixed by a 1 p, so all cycles of a 1 ^ p have length 2d. It follows 
that Bij = mn/(2d) = mgcd(j,n)/2. In case (m), two rows are fixed by a l p, 
so cycles of a 1 ^ p in these rows have length d and all others have length 2d. We 
conclude that B^ = 2n/d + (m — 2)n/ (2d) = (m + 2) gcd(j, n)/2. If the formula 
for B^ that holds when d is even were valid generally, we would have the second 
term in the formula for 62(771, n). The third term in the formula for bi(m, n) is 
a correction to the second term to treat the cases (i)~(iii) in which d is odd. 

Next, the formula for 63(771,71) follows by symmetry. More explicitly, 

m-ln-l n — 1 m — 1 

"»(•». ») = sb £ E 2C " = £ E 2 ""= m», ™> 

i=0 3=0 j'=0 i=0 

because the number of cycles Cij of c^t^B for an m x n array is equal to the 
number of cycles Bji of a J T l p for an n x m array. 

Finally, we consider 64(777, 77). For i = 0, 1, . . . , m — 1 and j = 0, 1, . . . , n — 1, 
a l r J p8 has the effect of reversing the first 777 — i rows, reversing the last i rows, 
reversing the first n — j columns, and reversing the last j columns. If m and 77 
are odd, then there is one fixed point and (77777 — l)/2 transpositions, so Dij = 
(77777 + 1)/2 for all i and j, hence 64(777, 77) = 2( mn ~ 3 )/ 2 . If m is odd and 77 is even, 
then = 77777/2 for all i and even j and = mn/2 + 1 for all i and odd j. 
This leads to 64(777,77) = (l/8)[2 m ™/ 2 + 2 mn ' 2+1 ] =3-2 mn / 2 - 3 . If m is even and 
77 is odd, then = mn/2 for even i and all j and = mn/2 + 1 for odd i and 
all j. This leads to the same formula for 64(777, 77). Finally, if 777 and 77 are even, 
then Dij = mn/2 unless i and j are both odd, in which case D^ = mn/2 + 2. 
This implies that 64(777,77) = (l/4)[(3/4)2" m / 2 + (l/4)2 m,l / 2 + 2 ] = 7 • 2" m / 2 " 4 . 



7 



This completes the proof. 



□ 



It is easy to check that 6(1, n) reduces to (2). Table 4 provides numerical 
values of 6(m, n) for small m and n. 



Table 4: The number 6(m, n) of toroidal m x n binary arrays, allowing rotation 
and/or reflection of the rows and/or the columns, for m,n = 1,2, ... ,8. 



13 
36 
158 
708 
4236 



158 
1459 
14676 



708 
14676 
340880 
8999762 
245619576 



13 
237 
4236 
184854 
8999762 
478070832 
26185264801 



18 
687 
26412 
2445918 
245619576 
26185264801 
2872221202512 



30 
2299 
180070 
33888844 
6873769668 
1466114420489 
321686550498774 



30 2299 180070 33888811 6873760668 1166111120189 32 1686^0 198774 72057630729710704 
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